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- - - Abstract 

' We calculated the bolometric images of relativistic slim disks during radiation-pressure-driven thermal 

, instability. When the mass-accretion rate exceeds the critical one, an inner region of the standard accretion 

disk bursts to change to a slim disk state having a large scale height. That is, the inner region of the disk 
JH , becomes high temperature, and the thickness of the disk increases due to the increase of the radiation 

pressure. As a result, we found that the observed image of the disk during the burst strongly depends on 
the inclination angle. That is, radiation from the innermost disk would be occulted by the disk outer rim 
for high inclination angles {i ^ 70°). We also calculated the spectral energy distribution during thermal 
instability. The Wien peak of the spectrum of high inclination angles becomes softer than that of low 
inclination angles due to the geometrical thickness. From these facts, even if the burst happens in black 
^ ' hole candidates, we may not observe the burst when the inclination angle is large. We may suggest that 

, numerous luminous black hole candidates are still hidden in our Galaxy. 
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Introduction 



Observational features of microquasars are considerably similar to those of quasars. That is, both microquasars and 
quasars have relatively high luminosities, relativistic jets, and central black holes (Mirabel, Rodriguez 1999; Fender, 
Q ' Bclloni 2004 for recent reviews). Apart from the difference of black hole masses, an investigation of microquasars is 
^ ^ strongly linked to understandings of quasar activities. Hence, many astronomers are engaged in studying this field. 
c/2 ■ The galactic object GRS 1915-f 105 is one of the famous microquasars in our Galaxy owning to the unique obser- 
' vational behavior; i.e., the speed of jets is about 0.92 c, quasi-periodic, flare like bursts, and limit-cyclic luminosity 
variations (Mirabel, Rodriguez 1999; Belloni et al. 2000; Fender, Belloni 2004). Although curious characteristics in 
this object collect attention of many astronomers, the pattern of variability is too complicated to understand well 
(Belloni et al. 2000). The X-ray light curves of GRS 1915-1-105 in its high luminosity state show limit-cyclic bursting 
■ behavior on the timescale of a few tens of seconds. The typical timescale of the limit-cyclic phenomena corresponds to 
a viscous time scale in an accretion disk. Accordingly, the phenomenon is attributed to the radiation-prcssure-driven 
thermal instability around the disk inner region (Belloni et al. 1997a, 1997b). 

The theory of the radiation-pressure-driven thermal instability in accretion disks was investigated by many authors 
(Shibazaki, Hoshi 1975; Shakura, Sunyaev 1976). Honma et al. (1991) was the first to numerically demonstrate the 
limit-cycle oscillations between the two different states, based on a time-dependent calculation. Unfortunately, this 
kind of instability could occur only when the mass-accretion rate is larger in comparison with a critical one, ^ O.ILe, 
and a candidate where such an instability occurs was not found until microquasars were discovered. They show 
that an existence of an adequately radiation-pressure dominated region is essential for developing thermal instability. 
After the discovery of microquasar GRS 1915-1-105, Szuszkiewicz and Miller (1997, 1998) confirmed the results of 
Honma et al. (1991), and moreover, they showed the difference of viscosity parameters, a, during thermal instability. 
The emergent spectrum during thermal instability was studied by Zampiri et al. (2001). Watarai and Mincshige 
(2003a; hereafter WM03) explained the limit-cycle variations of the maximum temperature and the emitting radius in 
GRS 1915-1-105. Recentry, it is confirmed that 2D SPH calculations also reproduce the limit-cycle oscillation (Teresi 
et al. 2004a,b) When we compare an observation with a disk model, there remain two important elements to examine 
more carefully. They are "geometrical thickness of the disk" and "special/general relativistic effects". If the accretion 
rate increases, the thickness of the disk rapidly increases due to the radiation pressure within the disk (Shakura, 
Sunyaev 1973). In addition, the "inner edge" of the disk in the supercritical accretion decreases down to less than 3 
Tg (Watarai, Mincshige 2003b). Hence, wc should take into account special/general relativistic effects; i.e., transverse 
Doppler effect, gravitational redshift, and light bending via curvature of space-time. These two important elements 



2 



A. Kawata et al. 



are strongly related to a viewing angle from a disk to a distant observer. Dependence of inclination angles in different 
accretion rates was pointed out by Watarai et al. (2005), and the observed spectrum significantly changes with high 
inclination angles and high mass accretion rates. However, they do not consider in the time evolution of the disk. 

Accordingly, we tried to examine the observational spectra during thermal instability in WM03 with various incli- 
nation angles. Note that previous papers have never considered the geometrical and special/general relativistic effects 
simultaneously during thermal instability. Thus, we examine there two effects in this paper. 

In the next section, we briefly overview the mechanism of thermal instability and their characteristics using the 
results of WM03. We show the flux images and the X-ray spectra of the inner disk during thermal instability in section 
3. The comparison with the X-ray observation of the microquasar GRS 1915+105 is discussed in section 4. The final 
section contains our concluding remarks. 



2. Disk Structure during Thermal Instability 



2.1. Basic Equations 

The basic equations for the viscous accretion disk are shown in the textbook by, e.g., Kato et al. (1998). We here 
briefly summarize the basic equations. 

We use the cylindrical coordinates (r, z), where the z-axis is the symmetrical one of the disk. We adopt a 
pseudo-Newtonian potential, tp = —GM/{R — r^) with R = \/r'^ + as the effect of general relativity (Paczynsky & 
Wiita 1980); where is the Schwarzschild radius defined by — IGMjt? = 3 x 10^ (M/M©) cm. The disk is assumed 
to be axisymmetric. The effects of outfiow/jet or coronal dissipation are ignored in this calculation. 

We adopt one-zone approximation in the vertical (z) direction, and the physical quantities are integrated in the 
vertical direction; i.e., the surface density is S — J_^pdz = 21^ pH, and the height-integrated pressure is 11 = J_^pdz = 
2/N+ipi?. Here, p, p and H are the mass density, the total pressure, and the scale height, respectively. The coefficients 
In and /n-i-i were introduced by Hoshi (1977). In the vertical direction, the density and pressure are assumed to be 
related by the polytropic relation p cx p^+^Z^. We assign iV = 3 in all calculations {I3 = 16/35, = 128/315). The 
height-integrated equation of state is 

n = Hrad + Hgas = ^T^* ■ 2H + ^STc, (1) 

3 pmn 

where the first term on the right-hand side represents the radiation pressure (a and T^. are the radiation constant 
and the temperature on the equatorial plane, respectively) and the second term represents the gas pressure (fee is the 
Boltzmann constant, /2 = 0.5 is the mean molecular weight, and mn is the hydrogen mass). We assume that, during 
thermal instability, the hydrostatic equilibrium holds in the vertical direction of the disk, 

{2N + 3)^^H'nl, (2) 

where I^k denotes the Keplerian angular frequency under the pseudo-Newtonian potential. 

The mass conservation, the radial component of the momentum conservation, and the angular momentum conser- 
vation are respectively written as follow; 

|(,S) + |-(rI]..)=0, (3) 
|(rl]..) + §-ir^vl + .H) = - 4) + (l - ^) n, (4) 

^(r^Ez;^) + ^{r^j:vrv^ - r^T,^) = 0, (5) 

where Vr and v^p are the radial and azimuthal components of velocity, respectively, and related to the angular momenta 
by I = rv^ and Iyl = rvK = r^^K, respectively. 

The r(/5-component of viscous stress tensor in equation (5) is prescribed as 

Tr^= I = -a/3^H (6) 

J-H 

where /i is a parameter (0 ^ ^ 1) and l3 = Pg^s/p- In the limit of /i — > 0, we recover the ordinary a viscosity 
prescription (Shakura & Sunyacv 1973). If we take the limit of /i^ 1, the viscosity depends on the gas pressure only. 
Finally, the energy equation is 

^ f ^ ^ 



-(rEetot) + — [r(Eetot + n)i'r - rTr^v^] = -2rF (7) 
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in which advective cooUng, viscous heating, and radiative coohng are considered. The exphcit form of the total energy 
etot is 



Etot 



3(l-/3) + :^ + l 



^ + Uvf.+vl) + Mr) (8) 



where the first term on the right hand side is the internal energy of the gas (7 is the adiabatic index and we set 7 = 5/3 
in the present calculation) , the second term represents the kinetic energy, and the last term is the potential energy on 
the equatorial plane, -00 = ^GM / {r — r^). Radiative cooling flux per unit surface area in an optically thick medium 
is given by 

where R is the average opacity written by an optical depth r as follow, 

T = = (Kes + Kff )S (10) 

where Kos = 0.4 is the opacity of the electron scattering, ns = 0.64 x 10^^ pT~'^^'^ is the absorption opacity via thermal 
Bremsstrahlung, and p= (16/35)p and T = (2/3)Tc are vertically average density and temperature, respectively. 

We solve equation (3)-(7) by the modified Lax Wendroff method with artificial viscosity. The calculations are 
performed from the outer radius at 2000 Tg down to the inner radius at ^ 2.2rg. The total mesh number is 250. We 
fixed black hole mass to be m = 10 (m = M/Mq, where Mq is solar mass). The critical accretion rate. Merit, is defined 
by L^jf? , where is the Eddington luminosity and c is the speed of light. We define the dimensionless accretion 
rate to be m = AI / Mait = Mc^/Lb throughout the present study. 

2.2. What Is Thermal Instability? 

In this section, we briefiy review the disk structures during thermal instability. Generally speaking, thermal insta- 
bility in an accretion disk occurs in radiation-prcssure-dominated regions (Shibazaki, Hoshi 1975; Shakura & Sunyaev 
1976). 

In standard accretion disk regimes, where the mass-accretion rate is lower than the critical one, the viscous heating 
energy, Q^isi Urst transformed to thermal energy, and finally released to the energy as a form of radiation Qj!^^. 
That is, the accretion time of the disk gas is longer than the radiative cooling time, thus the disk keeps the equilibrium 
state (Shakura, Sunyaev 1973). 

However, when the mass-accretion rate at the disk inner region increases, the radiative cooling energy cannot balance 
with the viscous heating energy. This is because that the dependence of a mass-accretion rate of the radiative cooling 
is Qrad °^ M^^^, while Q^jg c>c M for the viscous heating. When the M exceeds the critical rate, no equilibrium solution 
exists (Kato et al. 1998). 

As the accretion rate increases, the radial velocity of the flow increases, and finally the emitted photons are trapped 
by the advective fiow. That is, the advective cooling begins to work, Q~^^ esc M^. After the mass falls into black hole, 
the disk cannot maintain the slim disk state, thus the slim disk transits to the standard disk again. This transition is 
known as a "limit-cycle oscillation" . During the limit-cycle oscillation, the scale height and temperature distribution 
of the disk dramatically varies. We adopt the data by WM03 with the parameter set /j, = 0.1, a = 0.1 and m = 10. Time 
evolution of the disk scale height and temperature distribution during the transition is shown in figure 1 (WM03). 
Then, we verify that the geometrical thickness rapidly increases around the disk inner region (< lOOrg) during the 
disk transition due to the appearance of the radiation-pressure-dominant region. The disk thickness is close to 45 
degrees (H/r ^ 1) at the peak. Therefore, we should pay special attention to the geometrical thickness of the disk. 



3. Observed Characteristics during Thermal Instability 



We anticipate that photons emitted from an inner region are blocked by the increased scale height. In addition, 
photons from the inner region are affected by the special/general relativistic effects due to strong gravitational fields. 
Therefore, we expect that the image of an accretion disk is strongly affected by the scale height. The geometrical 
effect was considered in Madau (1988) and Fukue (2000). However, the special/general relativistic effects were not 
considered in these. The geometrical and relativistic effects have been considered in Watarai et al. (2005). Therefore, 
we pay special attention to the effect of scale height and make the images of the disk during thermal instability. To 
confirm how we see the accretion disk, we adopt the " Ray- Tracing Method" (The details of this method is written in 
Luminet 1979; Fukue, Yokoyama 1988). 

Photons emitted from any points on the disk travel along the null geodesies to be received by the distant observer. 
After the ray arrives at the disk surface, we obtain the physical quantities, such as, the temperature, velocity field, and 
scale height, which are numerically-calculated disk data at the arriving point. Due to relativistic effects, the observed 
radiative flux Fobs and blackbody temperature Tobs related to the emitted radiative flux and temperature Tem as 
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Fig. 1. Disk evolution during thermal instability; the left and right panels represent the effective temperature profiles and the disk 
scale height, respectively. The sequences are before the thermal instability (thick solid line), at the peak (thick dashed line), just 
after (dotted line), 30 seconds after thermal instability. 

p 

using Lorcntz invariancc. Here, redshift factor z of the light emitted from a surface of the disk: 1 + z = iJcm/^-obs = 
L~^'yD~^ (see Kato et al. 1998), where L is the lapse function representing the gravitational redshift, 7 is the Lorenz 
factor denoting the transverse Doppler effect, and D is the Doppler factor expressing the longitudinal Doppler effect. 

3.1. Imaging of the Disk during Thermal Instability 

The images of bursting disks and black hole shadows during thermal instability arc presented in figure 2. These 
figures represent the time evolution of the disk for several inclination angles (i = 1°, i = 50°, and i = 80° from left to 
right), and of the states of the disk (before the burst to after the burst from top to bottom). 

Let us see the image for i = 1° (face on). As the mass-accretion rate increases, the effective temperature of the 
innermost region of the disk rises up to ~ 10''K, the high energy radiation is expected around the vicinity of the black 
hole. After the burst peak, the temperature in the disk innermost region descends rapidly. 

On the other hand, for an intermediate inclination angle (i = 50°), the emergent radiation from the left part, which 
approaches the observer with the relativistic speed near to the black hole is remarkably enhanced by the Doppler 
beaming effect (blue shift), while the right part is reduced (red shift). Moreover, the light ray originating from the 
opposite side of the disk is strongly bent by the gravity of the black hole before it reaches the observer. Therefore, the 
images of the accretion disk become asymmetric and distorted shapes. However, the effect of self-occultation is not 
remarkable in the case of z = 50°. The self-occultation effect remarkably appears in the case of a high inclination angle 
{i = 80°). The black hole shadow is completely hidden by the disk itself in the peak phase. This effect of geometrical 
thickness affects the luminosity and spectral energy distribution (SED) of the black hole system. We inspect this 
geometrical effect for the luminosity in the next subsection. 

3.2. Luminosity Evolution during Thermal Instability 

Figure 3 represents the bolometric luminosity during thermal instability. The luminosity strongly depends on the 
viewing angle due to the projection effect and self-occultation. When the viewing angle is i = 1°, only gravitational 
redshift dominates. Therefore, the luminosity of the disk varies from Lobs ~ 0.01-0.4 L^dd- This result is consistent 
with that of WM03. As the inclination angle increases, the Doppler beaming and transverse redshift begin to be at 
work. However, the reduction of observed area mostly affects and observed luminosity decreases. In the case of i ^ 50°, 
not only the effect of projection but also the self-occultation affect to the image of the disk. So the luminosity of the 
disk with i — 80° decreases much more than i = 50°. Moreover, due to the self-occulting, the profile of the light curve 
for the case oi i = 80° is different from and flatter than the other cases of low inclination angles. 

To investigate how much the luminosity is blocked by the self-occultation effect, we plot the luminosity as a function 
of inclination angles (figure 4) . In figure 4, we show the relation between the light curves of thermal instability at the 
peak and the influence of inclination angle. The case of cos i ^ 0.8 falls below the evaluation considering only the 
projection effect, i.e., cos i because the gravitational redshift dominates. In the case of 0.4 < cos i ^ 0.8, the disk 
luminosity exceeds the evaluation. This is because the Doppler beaming dominates compared with the gravitational 
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Fig. 2. Observed temperature distributions during thermal instability. From top to bottom, before thermal instability {t = 418 s), 
the peak of thermal instability {t = 432 s), immediately after thermal instability [t = 444 s), of 20 seconds after thermal instability 
(t = 455 s). Unit of the time-scale is second. From left to right; i = 1°, i = 50° and i = 80°. The number of meshes is 250 X 250, 
and the screen size is ISr^ X 15r„. Unit of X-Y axis is the Schwarzschild radius r„. 
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Fig. 3. Light curves during thermal instability for several inclination angles; j = 1° (solid line), i = 50° (dashed line), i = 80° (dotted 
line). The abscissa is time in units of seconds, and the ordinate is the bolometric luminosity normalized by the Eddington one. 
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Fig. 4. Relationship between the peak luminosity and inclination angles i. The abscissa is cosi, and the ordinate is the disk 
luminosity normalized by the Eddington one. The dotted straight line is the projection factor under purely the geometrical effect. 

redshift. However, in the case of cos i ^ 0.4, the geometrical effect of self-occultation dominants and the luminosity 
decreases. Here, we emphasize that it is important to consider the geometrical effect for high inclination angles. 

3.3. Spectral Evolution during Thermal Instability 

We also calculated the spectra of the accretion disk. The spectra during thermal instability are affected by 
Comptonization due to their high temperature and large optical depth for electron scattering. Here, we adopt the 
diluted blackbody approximation as the effect of Comptonization for simplicity. The flux is given by 

^''-74(lV^-^^(/^eff) (13) 

where ly, i?^, and / are the frequency, the Planck function and the spectral hardening factor (Ebisuzaki et al. 1984). 
The spectral hardening factor represents the ratio of the color temperature to the effective temperature. The actual 
value of / is assumed to / = 1.7 which obtained from the radiation transfer calculation in the vertical direction of the 
disk (Shimura, Takahara 1995). We note that our present equation (13) includes the relativistic effect, i.e., 1/(1 + z)'^. 

The calculated spectral energy distributions of the accretion disk during thermal instability are presented (in actual, 
we plotted = 2i^^^ • 2-Krdr) In figure 5, the spectrum shapes for i = 50° are similar to z = 1° case throughout the 
time evolution. Because the geometrical thickness does not contribute to the image of the disk so much, the difference 
does not appear in the spectrum. On the other hand, in the case of i = 80°, the self-occulted effect clearly appears in 
the SED. The high energy photons emitted from the innermost region arc occulted by the large geometrical thickness, 
and the amount of the photons which finally reaching an infinite observer are decreases, and the spectrum becomes 
soft. We therefore should pay attention to the inclination angle. 

In this paper, we ignore the contribution of the coronal component. If the high temperature corona expands in the 
vertical direction above the disk, large amount of photons emitted from the disk are scattered by the electrons. It 
is expected that the coronal component contributes to the high energy spectrum during thermal instability. Strictly 
speaking, we should solve the radiative transfer equation until the surface of the optical depth for electron scattering 
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Fig. 5. Spectral energy distributions for several inclination angles with different time; before thermal instability (thick solid line), 
at the peak of thermal instability (thick dashed line), after thermal instability (dotted line), 30 seconds after finishing thermal 
instability. The inclination angle are (a) j = 1°, (b) i = 50°, and (c) i = 80°. 

Tcs ~ 1 surface including the coronal structure. However, we ignore such a coronal structure in this study in order to 
avoid a confusion of the origin of the spectrum. 

4. Discussion 

Where Are Hidden Limit-cycle Busters? 

In the previous section, we have demonstrated that the observed flux from the inner region of luminous accretion 
disks strongly depends on the inclination angle. According to our results, if other black hole candidates have a high 
inclination angle, it is difficult to detect the significant light curve variations. Because the disk outer rim blocks the 
emission from the disk inner region during the burst. This statement has already been pointed out by Watarai et al. 
(2005). 

Recently, type 2 counterparts of narrow line Scyfcrt 1 galaxies (hereafter NLSls) have been discovered by several 
authors (Nagar et al. 2002; Dewangan, Griffiths 2005). Type 2 NLSls candidates, NGC5506, NGC7582, and NGC7314, 
show the signatures of rapid X-ray variability (a few x 10^ to 10* seconds) and the reprocessed emission lines. According 
to the unified model of AGNs, type 2 AGNs has high inclination angles (Blandford 1990). That is, our present 
calculations can apply for accretion disks in type 2 NLSls. As for type 1 NLSl, supercritical accretion scenario seems 
to be consistent with the observations because of their large luminosity even for small black hole mass (Mineshige 
et al. 2000; Collin 2002; Collin, Kawaguchi 2004). If we believe the supercritical accretion scenario, we can predict 
that even if type 2 NLSls are really undergoing the thermal instability, the observed luminosity docs not reach at the 
Eddington luminosity due to the sclf-occultation effect. Unfortunately, the decreasing of the luminosity is caused by 
a lot of reasons, i.e., absorption of interstellar matters or molecular torus moving around AGNs. Thus, it is difficult 
to discriminate the cause of the decreasing of the luminosity. However, if it is possible to combine another physical 
information, such as time evolution of the reprocessed lines, it provides us to constraint the origin of the absorbers or 
the geometry of the accretion disk. 

4-2. Observational Implications of Disk Line Spectrum 

Since the discovery of the broadened iron Ka line at 6.4 keV in Seyfert galaxy MGC-6-30-15, an iron Ka line is 
widely used as a probe of the special/general relativistic effects near the supermassive black hole (Tanaka et al. 1995). 
The origin of line emission is roughly explained by the reflection model, i.e., the reflection component within the 
accretion disk or above the disk atmosphere absorbs the radiation from the intense X-ray source, and reradiates X-ray 
fluorescence lines (Miniutti, Fabian 2004). Actually, observed line spectra in some AGNs are consistent with above 
scenario. 

While the spectrum of narrow line Seyfert 1 galaxy 1H0707-495 shows more peculiar behaviors, for example, the 
variability time scale is about factor of 4 shorter than the typical AGNs, the photon index is steeper than usual AGNs, 
and the deep edge feature at about 7 kcV (Boiler et al. 2002; Tanaka et al. 2004). Fabian et al. (2002) assumed 
multiple reflection sheets in the radiation-pressure dominated disk to explain the deep edge at about 7 keV in NLSl 
IH 0707-495. This object also shows a strong soft X-ray excess, and the observed disk temperature implies that the 
bolometric luminosity in 1H0707-495 is close to the Eddington limit (Tanaka et al. 2004). When the disk is shining 
at the Eddington luminosity, the state of the disk should be in the state of slim disk rather than standard disk. If so, 
the reprocessed gas may be brown away by the strong radiation from the accretion disk, finally it may be difficult to 
reproduce the observed line intensity. We therefore should take into account the disk model more seriously because 
the reprocessed emission is strongly influenced by the radiation field of the disk. 

The correlation between the continuum and line spectrum during thermal instability seems to be important for 
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investigations of the spectral features like the emission line or absorption edge. Because the emission from the disk 
inner region is blocked by the disk outer edge for high viewing angles. How much part is covered with the disk outer 
edge? How is the absorption edge made? How docs the line flux correlate to continuum flux? These questions remain 
as a question that we should be going to answer in the near future. 

5. Conclusion 

We calculated the bolometric flux images and the observed spectra during thermal instability with considering the 
geometrical and special/general relativistic effects. Our results indicate that the observed images strongly depend on 
the viewing angles. In the case of low inclination angle (i ^ 70°), the image of the disk does not change so much. Only 
Doppler boosting and special/general relativistic effects appear. So the distribution of luminosity and spectrum has 
the peak in the peak phase. 

In contrast, if the inclination angle of the system is high (i 70°), the observed luminosity during the burst does 
change by much due to the self-occultation. Therefore, these facts indicate that if these situations occur in the real 
objects some bursting black holes still lurk in our Galaxy. If the jet from the accretion disk in the high luminosity 
state is related to the relaxation oscillation of the disk, we could detect this bursting bevavior through the jet. The 
disk-jet interaction will be an intriguing issue in high luminosity objects (Fender, Belloni 2004) 

Moreover, if we predict the relationship between the continuum and the line during the disk evolution we may 
restrict the physical state of the object from the comparison with the observed data. The evolution of the line spectra 
is one of our next issues. 
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